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Abstract 



Accurate calculation of internal and surface scattering losses in fused silica microspheres is done. 
We show that in microspheres internal scattering is partly inhibited as compared to losses in the 
bulk material. We pay attention on the effect of frozen thermodynamical capillary waves on surface 
roughness. We calculate also the value of mode splitting due to backscattering and other effects of 
this backscattering. _ 

Optical microsphere resonatorsa working on whispering gallery modes combine several unique features - very high 
quality factor, small size and effective volume of field localization with low cost, that make them very attractive for 
' future applications in optoelectronics and measurement science. Such mkgoresonators may be for example flsed as 
q ■ interferometers and filters with record finesseatl, for QED experimentsHoO and for diode laser stabilization!]. The 
^ ' last application as for now looks the most intriguing as it allows to create centimeter size cheap tunable laser with 
, kilohcrtz linewidth. 

Microspheres are prepared from fused silica or other pure glassy medium by autoforming under the action of surface 
tension in the flame of a microburner or in the CO2 beam. 
("^ ■ The properties of the microspheres are analyzed theoretically and experimentally rather wide. The theory of 
Oh| whispering gallery modes is well known from dectrodynamicsEil and methods of optimal coupling with them were 
elaborated and experimentally confirmed!]. Intl nonlinear properties of fused silica microspheres were investigated 
theoretically and experimentally. Mechanisms limiting the quality factor of microspheres were outlined in several 
J> , papers. It was shown that the main factor preventing the obtaining and preservation of high Q-factor is surface 
IO ' atmospheric water adsorptionl! However, values as. .high as Q ~ 8 x 10 9 - very close to fundamental limit of internal 
, losses were obtained at He-Ne and near IR regiorou. Attempts to obtain higher Q moving to minimum of losses of 
■ fused silica 1.55/xm, were as for now unsuccessful. Besides chemosorbed water the reason may be surface scattering 
on inhomogeneities. 

Problems of scattering inside the microsphere and on its surface were analyzed formerly incompletely. Scattering 
I leads not only to the limitation of the Q-factotliiit also lifts degeneracy between degenerate sine and cosine modes. 
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This effect may be observed as mode doubletsErlij. Especially badly explored are questions of surface scattering. In 
' different papers one may find expressions leading not only to different numerical estimates but to different functional 
dependencies from the size of the resonator and wavelengthlj 

v 

Oh 1- Scattering on internal thermodynamical inhomogeneities and quality-factor of microspheres 



Intrinsic scattering and absorption losses in microresonators were estimated previously from the bulk losses as: 



& : Q - ^ (i) 

where n is index of refraction a is intensity attenuation coefficient and A is wavelength. However, this approach is not 
quite accurate for scattering. We remind here the method of derivation of scattering coefficient a (see for exampleEj) 
to see what modifications should be made in Eq. (Q) to take into account specific features of the microspheres. 

Let us divide the whole volume of medium in small volume dv, each having due to fluctuations dielectric constant 
e(r) = 5e(f) + e . Internal small inhomogeneities in the field of the mode behave as dipoles reradiating light in all 
direction according to the Rayleigh formula: 



I s _ 7r 2 sin 2 tf 
7 ~ A 4 r 2 



5e(f 1 )<5e(f 2 )^ 1 ^ 2 ), (2) 
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where $ is the angle between the dipole axis (coinciding with polarization of field) and direction of scattering, r is the 
distance from scatterer. 

The next step is to integrate Eq.(0) over all angles on large sphere (r — > oo) to obtain total power of scattering. 



Ps = J J 5e(r 1 )5e(f 2 )dv 1 dv 2 ), (3) 

However, for the microsphere this is not correct. We should take into account total internal reflection (TIR). 

Beams falling on the surface under the angle larger than critical one 70 = arcsin(l/n) will either go back in the mode 
if this angle lies inside the mode's caustic or will be suppressed in destructive interference during several reflections. 
These beams may also go to another mode, which leads to internal mode coupling, but due to the rareness of mode 
spectrum this effect is negligible (specific case of coupling between oppositely circulating degenerate modes is analyzed 
below separately). In this way, only beams falling under angles less than critical should be added to losses. We may 
not take into account here Frensel transmission coefficients as these beams may leave the resonator during several 
reflections. Conditions for the cutting of angles for TE and TM modes are the following: 

/ T\ 2 / T\ 2 ^ 

sin 2 7 TS =(- ) (1 -sin 2 cos 2 <p) < — sin 2 j TM = ( ) sin 2 ?9 < — , (4) 
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here d is the distance of the dipole from the surface and a is the radius of the microsphere. If d <C a, that is always 
correct for high-Q whispering gallery modes, the first terms in (Q) may be omitted and hence the result will not 
depend from the size of the resonator. We skip here the following derivation of a by calculating thermodynamical 
calculations as conditions of angle cutting for large microspheres do not interfere with it (seeES). 

8ir 3 

otis = 7^n 8 p 2 nT/3 T , (5) 

where k is the Boltzman constant, T is the effective temperature of glassification (~ 1500K for fused silica), (3t is 
isothermic compressibility, and p is Pokkcls coefficient of optoclasticity at this temperature. 

Cutting conditions in this way may be taken into account by introducing suppression coefficients: 

2nn 

Qis = Kte,tm r (6) 

This coefficient Kte.tm is equal to the relation of complete scattered power to the power scattered on angles satisfying 
conditions (@). 

Numerical calculations for fused silica with n = 1.45 give 

Kte = 2.8 K TM = 9.6 (7) 

It follows from these values that TM-modes are less sensitive to intrinsic scattering losses, but these modes have 
stronger field on the surface, and therefore more sensitive to surface inhomogeneities and absorption on surface 
contaminations. 



2. Scattering on surface roughness 

To analyze surface scattering we analogously to the previous section calculate the value of a ss , describing losses of 
travelling wave per unit length. We start with the same expression in integral form (||), but now shall take into 
account only surface inhomogeneities. As before we should integrate this expression over angles with account of TIR, 
but for surface dipoles the part of light scattered above the surface may go free. In this way, suppression coefficient 
may be taken as 2K TE ^ TM /(K TE ^ TM + 1). In calculations of attenuation coefficient we as above for generality are 
not taking it into account and insert only in final formula for the quality-factor. 

Let the wave with intensity distribution I(y, z) travels along a guiding surface along the local x-axis, y-axis is 
chosen also along, and z-axis orthoganally to the surface. Small surface roughness leads to inhomogeneity of dielectric 
constant: 

Se(x,y,z) = (e - l)f(x,y)S(z), (8) 



2 



here 5(z) is delta function. If surface inhomogeneities are weakly correlated and their correlation function quickly 
abates to zero on the scale much smaller than the wavelength, roughness may be described by only two parameters - 
variance a = \J < f(x, y) 2 > and correlation length B. In this case 

/167T 2 
l(y,0)^ r {n 2 - l)nB 2 a 2 dx ^ Padx, (9) 

Thus, considering that the power of the wave is equal to P = J I(y, z)dydz, and considering that the wave travels 
close to the surface, we obtain: 

_ I(y,Q) 16(n 2 - 1)^ 2 ^ 



Jl(y,z)dz 3A 4 



Now we turn back to microsphere to calculate the ratio of intensities in above expression. For simplicity we limit 
ourselves only with TE-mode, which usually have higher quality factor. As the intensity is proportional to the square 
of electric field, and fields are described by Bessel functions, using the same approximations as irJ3 we obtain for 
TE-mode: 



dU(p)Y _ a(n 2 - 1) 



J ] 2 {knr)r 2 dr 

^ a 

a 2 ] 2 (kna) ~ 2] 2 (kna) \ dp ) p=kna ~ 2n 2 

Finally we obtain expression for the quality factor: 

V" - i + Kte 8 nn 2 B 2 a 2 { ' 

We may note that this expression is different from expression obtained in! but the reason is that,the authors under- 
estimated volumetric ratio [ll] considering it proportional to VaX and not a. In the same papem authors report on 
measurement of surface roughness of fused silica by means of scanning tunneling microscope. Values B = 5nm and 
a = 1.7nm were estimated. On Fig.l all cakylated limitations on quality factor are plotted, on this graphics UV and 
IR absorption in fused silica from literatureOj'EI are also taken into account. Bulk losses in fused silica were taken as: 

a ~ (0.7^m 4 /A 4 + 1.1 x 1(T 3 exp(4.6 fJ,m/\) + 4 x 10 12 exp(-56fim/\))dB/km (13) 
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FIGURES 
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Fig. 1. Influence of different types of losses on the quality factor of microspere 

For the estimates of surface scattering resonator radius a = was taken. One can see that for the measured 

Q in He-Ne range suppression of scattering is compensated by surface scattering. For longer wavelengthes, it seems, 
additional experiments are required. For very large spheres (several millimiters in diameter), quality factor sufficiently 
higher then 10 11 may be obtained— It seams that even Q ~ 10 12 is not impossible as Rayleigh scattering may be 
lowered by heat treatment on 25%E-X The only unknown factor is surface absorption on hemosorbed layers of OH~ 
ions and water. It should scale linearly with a as surface scattering (estimate inca is also based on the wrong volumetric 
ratio and in this way incorrect) but its dependence on A is unknown. It is known, however, that in the bulk hydroxil 
ions lead to vibronic absorption peak at A = 2.73/im and at obertones at 1.37, 0.95, 0.725 and 0.586^m. 

Nevertheless the problem of surface scattering in microsperes is not yet closed neither theoretically nor experimen- 
tally. It is reasonable to suggest that surface roughness may be attributed to surface capillary waves frozen during 
solidification. These waves leads to fluctuations: 

/(M) = 5>>l'(M), (14) 

where Y^f is spherical angular function, and L > 1. If according to thermodynamics each wave has energy kT, then0 
kT 

< hl >= *(L-l)(L + 2) ' (15) 

where a is the coefficient of surface tension ~ 200 dyn/cm for fused silica at temperaturer-jT = 1500-ftT. Though 
as estimates show the size of these fluctuations will be several times less than measured inB, correlation function, 
calculated for such inhomogeneities has logarithmic shape and in this way may not be characterized by correlation 
length. Fluctuations on the scale of wavelength are of the same order as in nanoscale. In other words, our approach will 



not work for this case. Unfortunately the estimate done ir£j fox-scattering on capillary waves Q cs ~ a\/((n — l)b L ) 
is also incorrect (it was shownt3 and experimentally confirmedtH that perturbation with M=0 (ellipticity if L=2), do 
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not perturb in the first order of approximation quality factor at all). The problem of frozen capillary waves in any 
case deserves special consideration. 



3. Coupled modes in microspheres 

Mode coupling in microspheres due to surface and internal inhomogeneities may be described using variational ap- 
proach. Random deviations of dielectric constant may be written in the form: 

Se = f(d,cj>)F(r), (16) 

where F(r) - is random radial function, and f(0, 4>) - is random angular function. In particular case of small surface 
roughness when random fluctuations of the surface of the sphere may be described as: 

r(M)=a + /(M), (17) 
and expression ( [l6|) may be written as 

5e = (n 2 - l)f(9^)S(r - a). (18) 

From the Maxwell equation, wave equation for the fields inside the microsphere with inhomogeneities may be 
obtained: 



The solutions of unperturbed equation without inhomogeneities (if Se — 0) have the form 
Ej = e ( - icjjt)e(r,6,<p), (20) 
where ej(r, (ft) - is vector harmonic, satisfying Helmholtz equation: 

Aej + e°k]ej = (21) 

and index j corresponds to all possible types of oscillations, and j = corresponds to initially excited mode. Using 
the method of slowly varying amplitudes we find solution as: 

E = e -^J2 A j(m (22) 
After substituting this sum in equation (|l9|) and omitting small terms we obtain: 



2w ° J2 + u 2 5eJ2Mm + ^°EK 2 - ofoMWi = (23) 

i 

After multiplication of this equation on ej* and integration over the whole volume, with account of modes' orthoganal- 
ity, we obtain usual equations for coupled modes: 

dAi- x — 

+ iAcu k A k = i Ajfrk, (24) 
j 

where Aujk = LUk — <-<Jo and 

Pjk ~ 2^ JWdv [ ' 

In this expression it is the random function Se which leads to the coupling between ej and e k - We are interested 
only on the modulus of the coefficient of 0jk, which determines the rate of energy redistribution between different 
modes. If the value of inhomogeneities and their correlation length are very small if compared with wavelength we 
may average f3 2 k and obtain that: 
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2 _ loI (J Se(r)Se(0)dv) 

^ ~ 4^ V~ k ' (26) 

where Vjk - is overlap volume of modes: 
_ /|ejf<fo/|efc| 2 (fo 

^" /|eJP|ei| a d« ■ ( } 

In the most interesting case of coupling between two modes A + (t) and A—(t), travelling inside the microsphere in 
opposite direction, fields' distribution for these two modes differs only on phase factor exp(±zm<^). In this way ej = 
and Vjk in this case transforms in effective volume of field localization: 

v *" ~ -JWj^> (28) 



4. Effects of coupled modes 



To analyze the consequences of internal coupling between modes travelling in opposite directions on the output 
characteristics of a resonator we may use the same quasigeometrical approach that we used inQ to analyze coupler 
devices for the microspheres. For simplicity we analyze here only the case of ideally mode matched (or monomode) 
coupler device. The set of equations for the internal and external amplitudes looks as follows: 

dv4- T 
—rr + i 6 o + S C + iAu)A + = iA_/3 + i—B m 

d^l_ 

—jT- + (5 + 5 C + iAuj)A- = iA + f3 



B t = y/l - T 2 B in + iTA+ 

B r = iTA- (29) 




Fig. 2. Backscattering in a microspere 

Where A + (t) and A- (t) are as before the amplitudes of oppositely circulating modes of TIR in the resonator (Fig. 
2) to model the whispering gallery modes. Bj„ is the amplitude of pump and Bt and B r are output amplitudes 
transmitted and reflected in coupler. T is the amplitude transmittance coefficient, describing coupler, 5q — 2ixnj a\ is 
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the decrement of internal losses, S c — T 2 /2tq is the decrement of coupler device, tq is the circulation time tq ~ 2irna/c, 
and Auj is frequency detuning from unperturbed resonance frequency ujo (for details seeo). The stationary solution 
of Eq. (p9h is the following 



A+ = 
A- = 

B t = 

B. r = 



T {S + S c ) 2 + [3 2 - Auj 2 + i2Auj{So + S c ) 
1 2S C {5 Q + S c + iAuj) 



Bj, 



T(S 
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- 5 C ) 2 +(3 2 - Auj 2 + i2Aaj(5 Q + 5, 
S 2 +(3 2 ~ Auj 2 + i2S Auj 



-Ba, 



Bj, 



{5a + ^ c ) 2 +(3 2 - Auj 2 + i2Auj{S + S e ) 
i2S c p 

~ {5o + 5 C ) 2 +[3 2 - Auj 2 + i2Auj{S a + 5 C ) 



B, 



(30) 



If internal mode coupling constant is weaker than attenuation (3 < Sq + 5 C than Eq.^ has only one resonance at 
Auj = and backscattering is small. The situation is not very different in this case from the case of one mode 
analyzed iro. In temporal language it means that internal coupling simply has no time to build backscattered wave 
during the ringdown time. Interesting is, however, that the regime of critical coupling (when B t — 0) is shifted and 
obtained not for 6 C = So but for S 2 = Sq + 1 and in this case not all input power is lost in the resonator but some 
part of it reflects back in the coupler. 
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(31) 



So + 5 C 

This means that ringdown time is equal to the time needed to repump circulating mode in oppositely circulating one. 




5 c /5 

Fig. 3. The dependence of power, reflected in coupler due to backscattering in microsphere, from loading 
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\B t ? 



(5c + S 



2 B in 



\B r 



{Sc + S ) 



:Bj. 



(32) 



The case of 8 > <5o + S c is much more interes ting and even le ads to somewhat unexpected results. In this case there 
are two resonances at frequencies Alu — ± V (0 2 — (8q + S c ) i.e. internal coupling lifts degeneracy between sine an 
cosine standing modes in the microsphere and there is enough time to form them. All intensities in resonances in this 
case do not depend on j3: 

What is interesting and inevident that if 5q *C S c (overcoupling) but still Sq + S c < 8 then most part of input 
power is backscattered and transmitted power tends to zero. This property may become extremely valuable for future 
applications of microspheres in laser stabilization. To verify this result additional experiments are required, but we 
saw many times that backscattering is practically absent when there is no splitting and practically does not depend 
on loading when doublets are clearly seen (see Fig. 3). 



5. Calculation of mode splitting on internal and surface inhomogeneities 



Mode coupling leads to splitting of initially degenerate modes, if 8 constant is much larger than mode decrement of 
internal and coupling attenuation <5o + S c then: 

Aw 28 

— = —, 33 

Ul LUq 

If thermodynamical inhomogeneities are calculated in the same way as before: 



Auj\ jn A p 2 KTp T I 3A 4 a 4S 



<■> J u V v *ff V 87rW *// 

in agreement with qualitative estimates ini. Effective volume of the most interesting T_E«i-mode may be calculated 
according to the formulaEl 

Veff = 2.3n' 7 / 6 a 11 / 6 X 7 / 6 , (35) 
In this way, for TiS«i-mode in fused silica microsphere: 
Aw\ 5 x 10" 7 /i™ 3/2 

~) ls ~ A n /i2 a 7/i2 ' ( 36 ) 
If I ^ m the following asymptotic approximation is valid: 

Veff, e m = V efm (l + 0.5V*- m- 0.5), (37) 

Now let us analyze the case of mode splitting due to the surface inhomogeneities. From (|25|) and (||) after averaging 
2 _ ul TTB 2 a 2 \e\ 4 

^ s " An*V eff J\e(r)\*dr (38) 



or for T^^!-mode in fused silica microsphere: 
Auj\ l.laB 



W J.. AV4 a 7/4' 



(39) 



It is easy to show that for measured size of surface inhomogeneitiesi this expressions gives substantially lower level 
of coupling between modes than internal inhomogeneities (36). On (Fig. 3) results of calculation for fused silica 
microsphere for A = 0. 63 fim according to (p6j,p9|) are shown. 
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Fig. 4. Normal mode splitting on inhomogeneities in microsperes 
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